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Abstract 

The current status of the determination of corrections to the hyperfine splitting of the ground 
state in hydrogen is considered. Improved calculations are provided taking into account the most 
recent value for the proton charge radius. Comparing experimental data with predictions for the 
hyperfine splitting, the Zemach radius of the proton is deduced to be 1.045(16) fm. Employing 
exponential parametrizations for the electromagnetic form factors we determine the magnetic radius 
of the proton to be 0.778(29) fm. Both values are compared with the corresponding ones derived 
from the data obtained in electron-proton scattering experiments and the data extracted from a 
rescaled difference between the hyperfine splittings in hydrogen and muonium. 
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I. INTRODUCTION 



High-precision measurements and calculations of energy spectra of hydrogen-like atoms 
provide tests of quantum electrodynamics (QED) with very high precision (see [l[Q, HQ] and 
references therein). In some cases, the current accuracy of QED calculations exceeds those 
of the known values of fundamental physical constants. For instance, recent measurements 
and calculations of the g factor of hydrogen-like carbon and oxygen have provided the basis 
for a new determination of the electron mass (see [5( and references therein). Measurements 
of the Lamb shift in hydrogen, combined with corresponding calculations, have facilitated 
to determine the Rydberg constant and to deduce an improved value for the proton charge 
radius flflE, <J. 

The relative experimental accuracy of the ground-state hyperfine splitting in hydrogen is 
better than 10" 12 mi. The error associated with the QED corrections to the hyperfine split- 
ting is estimated to contribute on the level 10 -9 . The major theoretical uncertainty arises 
from nuclear structure-dependent contributions. The most important structure-dependent 
term is the proton-size correction, which is determined exclusively by the spatial distribu- 
tions of the charge and the magnetic moment of the proton. It contributes on the relative 
level 10 -5 . Assuming that all other theoretical corrections are accurately known, one can 
determine the proton-size contribution by comparing theoretical and experimental values 
for the hyperfine splitting in hydrogen. The major goal of the present paper is to determine 
the Zemach and the magnetic radius of the proton by such a comparison. 

In Sec. II, we consider various theoretical contributions to the hyperfine splitting and 
derive the proton-size correction comparing theory and experiment. In Sec. Ill, we refine 
the value of the proton-size correction by recalculating some of these contributions and 
determine the magnetic radius of the proton, employing an exponential parametrization 
for the electric and magnetic form factors. The recalculation of the recoil correction has 
_d the va,ue of the Zemach radius conned to the previous result obtarned in 0. 
In Sec. IV, the results obtained are compared with corresponding data derived from elastic 
electron-proton scattering experiments 12, 13, Q| and the data extracted from a rescaled 
difference between the hyperfine splittings in hydrogen and muonium jis| . 

The relativistic units h = c = 1 are used throughout the paper. 



2 



II. HYPERFINE SPLITTING IN HYDROGEN 



The hyperfine splitting of the ground state in hydrogen can be written in the form 



AR 



theor 



E F (1 + <5 Dirac + 5 QED + 5 : 



structure > 



(1) 



where E F is the Fermi energy 16] 

g 

E F = -a(aZ) 3 



fh 

(m e + m p ) 3 /i N ' 
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/i p is the magnetic dipole moment of the proton, /xn is the nuclear magneton, m e and m p 
are the electron and proton mass, respectively. The relativistic correction <5 Dirac can easily 
be obtained from the Dirac equation [17]: 
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Here and in what follows we keep the nuclear charge number Z to separate the relativistic 



and radiative corrections 
^qed we re f er £ R e f s . 
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or recent achievements in calculations of the radiative correction 
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231 ] . The uncertainty of S^ ED is mainly determined 



by uncalculated terms of order a 3 (aZ) and by uncertainties associated with some of the 
calculated terms. The structure-dependent correction <5 structure i s usually expressed as the 



sum 
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The part associated with intrinsic proton dynamics (polarizability) 5 po1 has been recently 



evaluated in 



24j employing experimental and theoretical results for the structure functions 



of polarized protons. The correction due to muonic vacuum-polarization 5 Mvp has been 
obtained in 125], while the hadronic vacuum-polarization contribution 5 hvp was evaluated in 
HQ. ForUie weak interaction term 5 weak we refer to Refs. QQ. Va.ne 8 for all tne S e 



corrections together with corresponding uncertainties are presented in Table I. The Fermi 
energy E F is evaluated employing the values of the fundamental constants tabulated in 
The leading chiral logarithms contributions to the structure-dependent correction have been 
also investigated within an effective field theory j^. 

Now let us turn to the term S ngld , which is determined by electric and magnetic form 
factors of the proton. This quantity can be decomposed into two parts: <5 ngld = (5 ps + 5 reco11 , 



where 6 ps represents the proton-size correction and <5 reco11 is associated with recoil effects. The 
recoil part contains both terms arised from a pointlike Dirac proton and additional recoil 
correction due to the internal proton structure. Following Ref. [jy] we do not separate 
them. Calculations of the dominant contribution (relative order (aZ)m c /m p ) to the recoil 



correction have a long history (see |31[ and references therein). The contribution of the 
order (aZ) 2 m e /m p has been first derived in 13111. while the radiative-recoil correction of 
the order a(aZ)m e /m p has been obtained in [25J . To determine the magnetic radius of 
the proton from the hydrogen hyperfine splitting we propose the following. At first we 
calculate the structure-dependent part of the recoil correction in a rough approach, taking 
the proton magnetic radius to be the same as the charge one. Then we find the proton-size 
correction from a comparison of the experimental and theoretical values of the hyperfine 
splitting. Using the dipole parametrizations of the form factors we extract a preliminary 
value for the proton magnetic radius. Then we recalculate the recoil-structure correction 
with the obtained value of the proton magnetic radius and take into account the radiative 
and binding contributions to the proton-size term. At last we again find the proton-size 
correction and extract the magnetic radius. 

At first iteration we have calculated the recoil-structure correction (integrals VO, VV, k 2 , 
No.l, and No. 2 of Ref. [sjl) with the new proton charge radius (r 2 )^ 2 = 0.8750(68) fm [f| 
and with the same value for the magnetic radius of the proton. The total recoil correction 
turns out to be 5.84 ppm. Accordingly, the proton-size correction 5 ps is given by 

S ps = AE exp /E F - 1 - 5 Dirac - S QED - <f ccoil - S po1 - 5^ p - 5 hvp - <T eak . (5) 

This equation yields the value 5 ps = —39.98(61) x 10~ 6 . The next section is devoted to the 
deduction of the magnetic radius of the proton from the proton-size correction and to the 
next iteraction. 

III. MAGNETIC RADIUS OF THE PROTON 

The proton-size correction has been first evaluated in the non-relativistic limit by Zemach 

(5 Zemach = -2aZ memp R p} (6) 
m e + m p 
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where 

R p = [ d 3 r dV p E (r)p M (r') |r - r'| = -- f°° ^ G E (Q 2 ) G M (Q 2 ) - ll (7) 

J 71 Jo Q Z L /i p J 

defines the Zemach radius of the proton. Here Pe(f) and PmM denote the nuclear charge 
and magnetization distribution, respectively, where both densities are normalized to unity. 
The quantities Ge/m represent the electric and magnetic form factors, respectively. Since 
we deal here only with the static limit (Q° = 0), we define them to be dependent only on 
the spatial momentum transfer (squared), Q 2 > 0. The charge and magnetic mean-square 
radii are defined by the formula 

(r 2 ) E /M = / d 3 rr 2 p E/M (r). (8) 

In a first step we can approximate the proton-size correction by the Zemach formula. Then 
one can easily find the Zemach radius R p = 1.058 fm. The usual experimental fit for the 
proton form factors is the dipole parametrization 

which comes from the exponential model of the charge / magnetization distribution. But 
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35 



that for Q > 1 GeV 



recently, it was obtained in Jefferson Laboratory (JLab) 
the behavior of the electric form factor differs from the dipole parametrization. However, 
the Zemach correction is not sensitive to the form factors behavior for Q > 0.8 GeV. As it 
was shown in Ref. the contribution to the Zemach correction from the region Q > 0.8 
GeV is the same for different experimental models of the proton electric and magnetic form 
factors. Therefore, in what follows, we use the dipole parametrizations for the form factors. 
The parameters Rb, one for the electric and another for the magnetic form factor, are 
directly connected with the corresponding values of the root-mean-square radii. For the 
charge radius of the proton we take the value (r 2 )^ 2 = 0.8750(68) fm p obtained from the 
latest comparison of the theoretical and experimental values of the Lamb shift in hydrogen. 
Fixing the charge radius, we fit the magnetic radius such as to reproduce the Zemach radius. 
As the result we find a preliminary magnetic radius of the proton: (r 2 )^ 2 = 0.800 fm. In 
order to estimate the error associated with the model dependence we consider also the model 

G JL ab(Q 2 ) = (l " O- 13 ^) G d(Q 2 ) , (10) 



which is known as JLab model |34|. The error appeared is on the level of about 0.75%. 
In a second step we account for corrections to the Zemach formula 

^ps ^ Zemach _|_ ^radiative _|_ ^relativistic (11) 

where <5 radiatlve i s the radiative structure-dependent correction obtained in [2| and <5 relatlvlstlc 
is the binding correction derived in j3JJ. The radiative correction has been derived assuming 
the exponential model with the same parameter Rd for both charge and magnetization 
distributions, i.e. 



^radiative ^Zemach ^ 
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The accuracy of this approximation is sufficient for our purpose. Calculating <5 radiatlve for 
different i?o, we obtain <5 radiatlve = 0.0153(2) x <5 Zcmach . The binding correction has been 
expressed in terms of electric and magnetic moments of the proton: 



rrelativistic rZcmach 
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where 7 is Euler's constant and (r n )E/M = J d 3 r r n pE/M(^)- In part, this equation has 
been derived for the homogeneously charged sphere model for the proton charge distri- 
bution (with Rq = y / 5/3(r 2 )J/ 2 ). Nevertheless, the error induced by using this formula 
in comparison with other models for the charge distribution does not exceed 5%. Em- 
ploying the exponential model for the charge and magnetization distributions with electric 
and magnetic radii, (r 2 )^ 2 = 0.875 fm and (r 2 )]^ 2 = 0.800 fm, respectively, we obtain 
Rrelativistic = q.0002 x £ Zemach + 1.4 x 10" 8 . Thus the proton-size correction takes the form 

5 ps = 1.0154(2) x <5 Zcmach + 1.4 x 10~ 8 . (14) 

In addition, we need to correct the dominant term of the recoil contribution with the mag- 
netic radius (r 2 )^ 2 = 0.800 fm. As a result, the recoil correction turns out to be 5.94(6) 
ppm. 

Deducing again the proton-size correction by means of equation (j^J), we find the Zemach 
radius of the proton: 

i? p = 1.045(16)fm. (15) 



This value differs from the one obtained in 11], R p = 1.037(16) fm, mainly due to the 
recalculated recoil corrections with the new more precise charge and magnetic moment 
distributions. 



IV. DISCUSSION 

The value for the Zemach radius obtained above enables us to determine an improved 
magnetic radius of the proton: 

(r 2 )^ /2 = 0.778(29)fm. (16) 

The corresponding uncertainty is mainly due to errors associated with the polarizability 
effect as well as the uncertainty of the charge radius of the proton. In Table I we present 
the final values for the contributions to the hyperfine splitting in hydrogen. The value for 
S ps has been obtained by means of the experimental energy splitting, according to equation 

Another way to determine the proton magnetic radius is based on experimental data from 
elastic electron-proton scattering. Accordingly, Friar and Sick have recently determined the 
Zemach radius of the proton to be R p = 1.086(12) fm [3] and the proton charge radius 
(r 2 )j/ 2 = 0.895(18) fm |l2j. Based on these values for R p and (r 2 )^ 2 and employing an 
exponential parametrization for both electric and magnetic form factors, we find the value 
of the proton magnetic radius (r 2 )^ 2 = 0.824(27) fm. This value is in a good agreement 



with the recent experimental value of Sick - 0.855(35) fm presented in 14J, and also the 
result of Hammer and Meifiner 0.857 fm 141 is not far away. 

Recently Brodsky et al. [15J proposed another method to extract the Zemach radius. 
They considered a rescaled difference between the hyperfine splittings in hydrogen and muo- 
nium 

a^w^ = 1 + 5 ' (17) 

where p and /i indicate quantities which refer to the proton and muon, respectively, AE£ xp = 



4 463.302 765(53) MHz 38j, and Ep is the Fermi energy for muonium. Employing recent 
values of the fundamental constants |(| they have obtained 5 his = 145.51 ppm Q]. The 
ground state hyperfine splitting in muonium can be written as 

AE? heoT = E§(1 + 5 Dirac + 6 QED + Sl econ + 5j vp + 5" cak ) . (18) 



The corrections <5 Dirac and <5^ ED are the same as in the case of hydrogen, £^ vp and <5™ eak are 
the hadronic vacuum-polarization and the weak interaction contributions, respectively, and 
^recoil j g reco j} term, which consists of relativistic and radiative parts. From this formula 
together with Eqs. ([TJ and (|17|) one can immediately derive the proton structure correction 

^structure -)- ^ recon _|_ ^hvp _|_ ^weak _|_ ^hfs/^Dirac _|_ ^QED _|_ ^recoil _|_ ^hvp _i_ ^weak"\ (19) 

'. f, followin g to Ref . , we take into account only the relativistic part of the recoil correction 
£l| and neglect the contributions 5j vp and 5T ak , we obtain «^ tructurc = -31.8 
ppm. This yields the Zemach radius R p = 1.019(16) fm |l5( that differs significantly from 
our result, R p = 1.045(16) fm. We have found, however, that the difference disappears 
if one includes the omitted terms. This is mainly due to the radiative-recoil correction 



41 



42 



43 



With this term included, the total recoil correction 



evaluated in |20j, |3JJ |4C 

is determined a s # re co11 = —178.33 ppm. The hadronic vacuum-polarization contribution 



obtained in j3, I^J is £^ vp = 0.05 p pm , while the value of the correction due to Z°-boson 



exchange yields <5™ eak = —0.01 ppm [28, 



•-<■•-<■■- r: - - i- j, |29|. Substituting these values into expression (|T9"j) . 
we find ^tructure — —32.64 ppm. Utilizing the values presented in Table I, we obtain for the 
proton-size correction 5 ps = —40.15 ppm. Then the Zemach radius can be easily determined 
with the result R p = 1.047(16) fm, which is very close to the value obtained in this work, 
R p = 1.045(16) fm. 

As one can see, only a disagreement with the value for the Zemach radius and, therefore, 
with this for the magnetic radius, as is obtained from the electron-proton scattering experi- 
ments remains. At present we have no explanation for this deviation. One may hope, that a 

new determination of the proton charge radius via the Lamb shift experiment with muonic 

"1 

hydrogen, which is now in progress at PSI (Paul Scherrer Institute) |45j, will elucidate the 
situation. From the theoretical point of view, an independent calculation of the proton 
polarizability effect would be also desirable. 
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TABLE I: Numerical values for various corrections to the hyperfine splitting in hydrogen together 
with the assigned errors. The energies AE exp and Ep are given in units of MHz. 





Value 


Error 


Ref. 




1 420.405 751 767 


0.000 000 001 [10J 


Ep 


1 418.840 08 


n nnn 09 


ID) 


AE exp /Ep 


1.001 103 49 


n nnn nnn m 

u.uuu uuu u± 




^Dirac 


0.000 079 88 




[1 71 


^QED 


0.001 056 21 


0.000 000 001 


[1£, 12, 20, 21, 22,23J 


£PS 


- 0.000 040 11 


0.000 000 61 




^recoil 


0.000 005 97 


0.000 000 06 


[25. 311 , this work 


S pol 


0.000 001 4 


0.000 000 6 


[24] 


S^P 


0.000 000 07 


0.000 000 02 


[25] 


(jhvp 


0.000 000 01 




[26, 27] 


^weak 


0.000 000 06 




[2a, 29J 
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